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Abstract 

One of remarkable features of the gamma ray blazar Markarian 501 is the reported 
shape of the TeV spectrum, which during strong flares of the source remains essen- 
tiahy stable despite dramatic variations of the absolute 7-ray flux. I argue that this 
unusual behavior of the source could be explained assuming that the TeV emission 
is a result of synchrotron radiation of extremely high energy {E > 10^9 eV) protons 
in highly magnetized (5 ~ 30 — 100 G) compact regions of the jet with typical size 
R ^ 10^^ — 10^^ cm and Doppler factor ~ 10 — 30. It is shown that if protons are 
accelerated at the maximum possible rate, i.e. tacc = vii^g/c) with so-called gyro- 
factor r/ ~ 1, the synchrotron cooling of protons could not only dominate over other 
radiative and non-radiative losses, but could also provide good fits (within uncertain- 
ties introduced by extragalactic 7-ray extinction) to the 7-radiation of two firmly 
established TeV blazars - Markarian 501 and Markarian 421. Remarkably, if the 
proton acceleration takes place in the regime dominated by synchrotron losses, the 
spectral shape of the Doppler-boosted 7-radiation in the observer's frame is deter- 
mined essentially by the self-regulated "synchrotron cutoff" at eo — 0.3 6^r]~^ TeV. 
The hypothesis of the proton-synchrotron origin of TeV flares of BL Lac objects 
inevitably implies that the energy contained in the form of magnetic field in the 
7-ray emitting region exceeds the kinetic energy of accelerated protons. 
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1 Introduction 



Blazars are Active Galactic Nuclei (AGN) dominated by a highly variable com- 
ponent of non-thermal radiation produced in relativistic jets close to the line 
of sight (e.g. Begelman et al. 1984, Urry & Padovani 1995). The dramatically 
enhanced fluxes of the Dopplcr-boostcd radiation, coupled with the fortuitous 
orientation of the jets towards the observer, make these objects ideal labora- 
tories to reveal the underlying physics of AGN jets through multi-wavelength 
studies of temporal and spectral characteristics of radiation from radio to very 
high energy 7-rays (Ulrich et al. 1997) . First of all this concerns the BL Lac- 
ertac (BL Lac) objects - a sub-population of blazars of which two prominent 
representatives, Markarian 421 and Markarian 501, are firmly established as 
TeV 7-ray emitters. The flux variability on different time-scales (and, plau- 
sibly, of different origin) is a remarkable feature of TeV radiation of BL Lac 
objects. It ranges from the spectacular 1996 May 15 flare of Markarian 421 
with duration less than 1 h to the extraordinary high state of Markarian 501 
in 1997 lasting several months (for review see Aharonian 1999; Catanese & 
Weekes 1999). The recent multi- wavelength campaigns revealed that the TeV 
flares of both objects, Markarian 501 (Plan et al. 1998; Catanese et al. 1997; 
Krawczynski et al. 2000; Sambruna et al. 2000) and Markarian 421 (Buckley 
et al. 1996; Maraschi et al. 1999; Takahashi et al. 1999), correlate with X- 
radiation on time-scales of hours or less. This is often interpreted as a strong 
argument in favor of the so-called synchrotron- Compton jet emission models 
in which the same population of ultra-relativistic electrons is responsible for 
production of both X-rays and TeV 7-rays via synchrotron radiation and in- 
verse Compton scattering, respectively (see, e.g., Ulrich et al. 1997). However, 
the very fact of correlation between X-ray and TeV 7-ray fluxes docs not yet 
rule out other possibilities, in particular the so-called hadronic models which 
assume that the observed 7-ray emission is initiated by accelerated protons 
interacting with ambient gas or low-frequency radiation. 

Generally, the hadronic models do not offer efficient 7-ray production mecha- 
nisms in the jet. For example, for any reasonable acceleration power of protons, 
Lp < lO^^erg/s, the density of the thermal plasma in the jet should exceed 
10^ cm"^ in order to interpret the reported TeV flares of Markarian 501 by 
7r°-decay 7-rays produced at p-p interactions. Therefore this mechanism could 
be effectively realized only in a scenario like "relativistic jet meets target" 
(Morrison et al. 1984), i.e. assuming that 7-radiation is produced in dense 
gas clouds that move across the jet (e.g. Dar & Laor 1997) Recently a novel, 
"non-acceleration" version of 7r°-decay 7-ray production by blazar jets was 
suggested by Pohl & Schhckeiser (2000). 

The Proton Induced Cascade (PIC) model (Mannheim 1993; Mannheim 1996) 
is another attractive possibility for production of high energy 7-rays. This 



2 



model relates the observed 7-radiation to the development of pair cascades 
in the jet triggered by secondary "photo-meson" products (7-rays and elec- 
trons) produced at interactions of accelerated protons with low-frequency syn- 
chrotron radiation. The efficiency of this model significantly increases with 
energy of accelerated protons, therefore the postulation of an existence of ex- 
tremely high energy (EHE; E > 10^^ eV) protons is a key assumption for the 
PIC model. In a compact 7 production region of the jet with characteristic size 
less than 10^^ cm, the protons could be accelerated to such high energies only 
in the presence of large magnetic field, B IG. Even so, below I will show 
that the very fact of observations of multi-TeV 7-rays from Markarian 421 and 
Markarian 501 allow a rather robust lower limit on the "photo-meson" cool- 
ing time of protons, tp^ > 10^ s. This implies uncomfortably high luminosity 
in EHE protons which would be required to match the observed TeV 7-ray 
fluxes. 

Meanwhile, at such conditions the synchrotron radiation of the EHE protons 
becomes a very effective channel of production of high energy 7-rays. In this 
paper 1 show that for a reasonable set of parameters, which characterize the 
small-scale (sub-par sec) jets in Markarian 421 and Markarian 501, the syn- 
chrotron radiation of EHE protons not only may dominate over other possible 
radiative and non-radiative losses, but also could provide adequate fits to the 
observed TeV spectra of both objects^ . Moreover, this hypothesis could natu- 
rally explain one of the remarkable features of TeV flares of Markarian 501 - its 
essentially stable spectral shape despite spectacular variations of the absolute 
TeV flux up to factor of 10 or more on time-scales less than 1 day. 



2 Synchrotron radiation of protons 

The comprehensively developed theory of synchrotron radiation of relativistic 
electrons (see, e.g., Ginzburg & Syrovatsky 1965; Blumenthal & Gould 1970) 
can be readily applied to the proton-synchrotron radiation by re-scaling the 
Larmor frequency z/l = eB/2nmc by the factor 171^/171^ — 1836. For the 
same energy of electrons and protons, Ep = = E, the energy loss rate 
of protons {dE/dt)sy appears (mp/me)^ — 10-*^^ times slower than the energy 
loss rate of electrons. Also, the characteristic frequency of the synchrotron 
radiation z/c = 3/2 i/l {E/mc^)^ emitted by a proton is (mp/me)^ — 6-10^ times 
smaller than the characteristic frequency of synchrotron photons emitted by 
an electron of the same energy. Then the synchrotron cooling time of a proton, 

^ The synchrotron radiation of protons was actually included in the overall PIC 
code of Mannheim (1993), but the effect of this process was somehow disregarded. 
The importance of the proton synchrotron radiation in the PIC scenario was recently 
recognized, independent of the present paper, by Miicke &; Protheroe (2000). 
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tsy — E/{dE/dt)sy, and the characteristic energy of synchrotron photons ec = 
produced in the magnetic field B are 



and 

^ heBE 



where Sioo = -B/lOO G and £^19 = E/10^^ eV. Hereafter it is assumed that the 
magnetic field is distributed isotropically, i.e. B± — sin ip B with sin ip — y'^/S. 

The average energy of synchrotron photons produced by a particle of energy 
E is equal to em ~ 0.29 ec (Ginzburg & Syrovatsky 1965). Correspondingly, 
the characteristic time of radiation of a synchrotron photon of energy e by a 
proton in the magnetic field B is 

tsy(e) ~ 2.2 X 10^ B^oo^ (e/1 GeV)-^/^ ^ _ 



For comparison, the time needed for radiation of a synchrotron 7-ray photon 
of the same energy e and in the same magnetic field, but by an electron, is 
shorter by a factor (mp/me)^/^ — 1.5 x 10^. 

The spectral distribution of synchrotron radiation emitted by a proton of 
energy E is described by the equation 



where x — e/cc, and F{x) = x dxK^/^^x) ; K5/3{x) is the modified Bessel 
function of 5/3 order. The function F{x) could be presented in a simple ana- 
lytical form (e.g. Melrose 1980): 

F{x) = C x^/^ exp {-x) . (5) 



Numerical calculations show that with C ~ 1.85 this approximation provides 
very good, less than 1 per cent accuracy in the region of the maximum at 
X ~ 0.3, and still reasonable (less than several per cent) accuracy in the broad 
dynamical region 0.1 < x < 10. 
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Fig. 1. Different possible spectra of accelerated protons. At energies E <^ Eq all 
spectra have power-law behavior with Op = 2, but in the region of the cutoff Eq 
they have essentially different shapes (see the text). 

For the given energy spectrum Np{E) of relativistic protons, distributed isotrop- 
ically in a source at a distance d from the observer, the differential flux of 
synchrotron radiation is defined as 



—1 ^ 



(6) 



Figure 1 shows four examples of possible proton spectra. The curve 1 corre- 
sponds to the most "standard" assumption for the spectrum of accelerated 
(e.g. by shock waves) particles - power-law with an exponential cutoff at en- 
ergy Eq: 



(7) 



The curve 2 corresponds to a less realistic, truncated power-law spectrum, i.e. 
Np{E) oc at E < Eq, and Np{E) = Q at E > Eq. 

While the cutoff energy Eq in the spectrum of accelerated particles could be 
estimated quite confidently from the balance between the particle accelera- 
tion and the energy loss rates, the shape of the resulting spectrum in the 
cutoff region depends on several circumstances - the specific mechanisms of 
acceleration and energy dissipation, the diffusion coefficient, etc. For example. 
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the recently revived interest in the diffusive shock acceleration of electrons 
resulted in predictions that in the shock acceleration scheme one may expect 
not only spectral cutoffs, but perhaps also pronounced pilc-ups preceding the 
cutoffs (Melrose & Crouch 1997; Protheroe & Stanev 1999; Drury et al. 1999). 
Two such spectra are shown in Fig. 1. The curve 3 represents the extreme 
class of spectra containing a sharp (with an amplitude of factor of 10) spike at 
the very edge of the spectrum. The curve 4 corresponds to a smoother spec- 
trum with a modest pile- up (or "bump") and super-exponential (but not yet 
abrupt) cutoff. 




Fig. 2. The Spectral Energy Distributions of the synchrotron radiation produced 
by protons with spectra shown in Fig. 1. The curve 1 corresponds to the proton 
spectrum described by Eq.(7); the curve 2 corresponds to the truncated proton 
spectrum; the curve 3 corresponds to the proton spectrum with a sharp pile-up 
and an abrupt cutoff at Eq; the curve 4 corresponds to the proton spectrum with 
a smooth pile-up and a super-exponential cutoff. For comparison, the SED of the 
synchrotron radiation of mono-energetic protons, xF{x) oc x*^/^ exp (— x), is also 
shown (curve 5). 

The corresponding spectral energy distributions (SED: uS^ = e^J(e)) of syn- 
chrotron radiation are shown in Fig. 2. In the high energy range, e > eo, where 
eo is defined by Eq.(2) as eo = ec(£'o), the radiation spectrum from the proton 
distribution with sharp pile-up and abrupt cut-off is quite similar to the syn- 
chrotron spectrum from mono-energetic protons, xF(x) = x^^^ e~^; x = e/eo. 
This is explained by the radiation component associated with the line-type 
feature in the proton spectrum at E — Eq. 
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All synchrotron spectra shown in Fig. 2 exhibit, despite their essentially dif- 
ferent shapes, cutoffs approximately at x ~ f, if one defines the cutoff as 
as the energy at which the differential spectrum drops to 1/e of its extrapo- 
lated (from low energies) power-law value. Therefore the energy eo = €c{Eq) 
could be treated as an appropriate parameter representing the synchrotron 
cutoff for a quite broad class of proton distributions. In the case of mono- 
energetic protons, the cutoff energy coincides exactly with eo (see Eq. 5). This 
is true also for the power-law proton spectrum with exponential cutoff given 
by Eq. (7), for which the SED of the synchrotron radiation has a shape close 
to uS^ oc e^/^ exp [—{e/eoY^^] (see below). 

Below 1 will not specify the proton acceleration mechanisms, but rather as- 
sume, with some exceptions notified otherwise, that the accelerated protons 
in the source are represented by a "standard" featureless spectrum given by 
Eq. (7). For this energy distribution of protons, the delta- functional approx- 
imation gives a simple analytical expression for the differential spectrum of 
synchrotron radiation: J{E) oc exp [— (e/^eo)i/2]^ where F = (ftp + l)/2, 
and ^ is a free parameter introduced in order to adapt this formula to the 
accurate numerical calculations. 

The spectra of synchrotron radiation calculated in the delta-functional ap- 
proximation for q = 1/3, 2/3, 1, and4/3 (curves 2, 3, 4, and 5, respectively) 
are shown in Fig. 3, together with the result of accurate numerical calculations 
(curve 1). It is seen that the best broad-band fit, with accuracy better than 
25% in the entire region up to e ~ 20eo, is provided by = 1. Although — 2/3 
gives, in fact, better fit at e < eo, it significantly underestimates the flux above 
eo. And finally, both ^ = 1/3 and <; = 4/3 corresponding to the maximums in 
Si, and uSi, distributions of radiation emitted by a mono-energetic electron, 
do not adequately describe the radiation spectrum at e > O.leo. Thus, the 
most appropriate value for q is close to 1. This proves that the synchrotron 
cutoff indeed takes place at eo. At energies e <^ eo the radiation spectrum has a 
power-law shape with photon index («p + l)/2. Note that the steepening of the 
7-ray spectrum at higher energies is significantly smoother (oc exp [— (e/eo)^^^]) 
than the exponential cutoff in the spectrum of the parent protons ^ . For the 
power-law proton distribution truncated at Eq, the synchrotron spectrum be- 
yond eo extends approximately as exp [— e/eo]). 

Assuming now that the synchrotron radiation is emitted by a relativistically 
moving cloud of plasma (often called in literature as plasmon or blob), in which 
the spectrum of accelerated protons are represented by Eq. (7), the differential 
flux of radiation in the delta-functional approximation can be presented in the 



^ In fact, the numerical calculations show that the spectrum in the cutoff region 
drops even slower than exp [— (e/eo)^'^^]; compare the curves 1 and 4 in Fig. 3. 
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Fig. 3. Spectral Energy Distribution (SED) of the proton-synchrotron radiation cal- 
culated in the delta-functional approximation for the power-law spectrum of protons 
with a spectral index Op = 2 and an exponential cutoff at Eq (Eq. 7). The curves 2, 
3, 4, and 5 correspond to the parameter ? = 1/3,2/3, 1 and 4/3, respectively. The 
curve 1 correspond to accurate numerical calculations. 

following simple form 



where A = {Nq/AMH^) {a^B^ml/ml^) a("p"3)/25K+5)/2^ p _ (ap + l)/2. 



Co = ^j^cl-Eo), and 5j is the jet's Doppler factor. 

This convenient equation could be used for the fit of TeV observations from 
BL Lac objects in the framework of the proton-synchrotron model, assuming 
that the spectrum of the accelerated protons is given by Eq. (7) 



3 Efficiency of proton synciirotron radiation 

The efficiency of a 7-ray production mechanism is characterized by the ratio 

of the radiative cooling time to the typical dynamical time-scale of the source 
corresponding to the minimum radiation variability in the jet's frame 




(8) 



t* = At 



5j ~ 1.08 X 10^ Atsh^ios, 



(9) 



where Atsh — Ai/3h is the radiation variabihty in the observer's frame in 
units of 3 hour - a characteristic time-scale observed from Markarian 501 in 
TeV 7-rays (Aharonian et ah 1999a; Quinn et al. 1999), and 5iq = The 
proton-synchrotron radiation becomes an effective mechanism of 7-radiation 
with coohng time < 10^ s only when the photons are emitted by EHE protons 
with E > 10^^ eV in a strong magnetic field close to 100 G. In this regime, the 
synchrotron losses well dominate over non-radiative losses (see Fig. 4) caused 
by adiabatic expansion of the source or escape of particles from the source 
with characteristic time-scales (c/vq) x t* and (c/vesc) x t*, even in the case 
of relativistically expanding source {vq ~ c), or energy-independent escape of 
particles with speed of light, Vesc ~ c. 

It is interesting to compare the proton-synchrotron cooling time tgy with the 
photo-meson cooling time, tp^. The energy flux of low-frequency radiation that 
originates in the jet could be presented in the following form 

1/5^ = 10"^^ 5fir {hu/0.01 eV)-^+^ erg/cmh (10) 

where (yf^r is a scaling factor indicating the level of the far infrared (FIR) flux at 
hu = 0.01 eV (A ~ 100 /im) in unites of 10~^^ erg/cm^s. The flux normalization 
at FIR wavelengths, which play a major role in proton-photon interactions in 
the jet, makes the final results of calculations quite insensitive to the choice of 
the spectral index s. In the case of Markarian 501, the power-law extrapolation 
of the observed synchrotron flux at soft X-rays (< 1 keV) towards far infrared 
wavelengths with s = 0.5 gives (^fij. ~ 0.3. On the other hand, assuming that 
the whole radio-to-infrared flux observed from the direction of Markarian 501 
(see, e.g.. Plan et al. 1998) is produced in a small-scale (7-ray emitting) jet, 
the factor gf^j. could be one order of magnitude larger. 

For a source with the co-moving frame luminosity L\u')u, the observer at 
a distance d — cz/Hq would detect a flux (Lind & Blandford 1985) Si, — 
5^L'{v/5^)/At:(P (hereafter, the value Hq — 60km/sMpc will be used for the 
Hubble constant). The energy losses of protons in a low- frequency photon field 
with a broad power-law spectrum are dominated by photomeson-processes. 
For broad and flat spectra of target photons the photo-meson cooling time of 
protons can be estimated as tp^{E) ~ (c < cXp^/ > n(z/*)/iz/*)~^, thus for the 
flux given by Eq. (10) with s=0.5, 

V ^ 4.5 X 10^ A4 4^ {z/Om)-' g^l^ E^,'-' s , (11) 

where < (Jp^f >~ 10~^^cm^ is the photo-meson production cross section 
weighted by inelasticity at the photon energy ~ 300 MeV in the proton rest 
frame, and hiy* ~ 0.03Ef9^eV (see, e.g., Stecker 1968; Miicke et al. 1999). 
Thus, for both Markarian 501 and Markarian 421 with redshifts z — 0.034 and 
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z = 0.031, respectively, the photo- meson cooHng time cannot be significantly 
less than 10'' s, unless we assume a very high ambient photon density. Formally 
this could be possible, for example, adopting a much smaller blob size than it 
follows from the observed flux variability At^h ~ 1, and/or assuming a small 
Doppler factor of the jet, 5j <S 10. However, the photon density in the source 
cannot be arbitrarily increased, otherwise it would result in a catastrophic 
absorption of TeV radiation inside the source. 

In the field of ambient photons with differential power-law spectrum n{i') oc 
the optical depth for tlic photon-photon absorption is equal to Tyy{€) — 
A(s)(aT/2) huon{hi^o)R, where huo = A{m^c^y/e, and A{s) = 7/12 • 4^+^(5 + 
1)~^/^/ (s -|- 2) (Svensson 1987). For the spectral index s = 0.5, 

~ 0.4 Ai-^i 5ro' ^fir (z/Om)^ (e/1 TeV) V2 . (12) 



The optical depth Xy-y increases with energy as r^-y oc e^/^, therefore at 1 TeV 
it should not significantly exceed 1, otherwise the absorption of > 10 TeV 
7-rays becomes unacceptably large. For Markarian 501, assuming that the 
low-frequency radiation in the blob is described by Eq. (10) with s = 0.5 and 
adopting a rather relaxed estimate for g^^^ ~ 0.3, the condition r^^(l TeV) < 1 
results in a robust lower hmit on the blob's Doppler factor Sj > 7. 

10^ 
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Fig. 4. The characteristic acceleration and energy loss time-scales of protons in the 
blob. At each curve the scaling factors (the products of relevant physical parameters) 
are shown (see the text). 

It is worth noting that formally we may expect non-negligible flux of TeV 7- 
rays even at conditions when r-y-y S> 1, the radiation being contributed from the 
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"last layer" of the source with r^^ ~ 1. For a (quasi) homogeneous spherical 
source, this implies that the "TeV photosphere" occupies only part of 

the total volume of the source. On the other hand, the initial TeV radiation, 
produced throughout the entire source, is re-radiated in the form of secondary 
(cascade) electromagnetic radiation in the energy band below e*, where e* is 
determined by the condition r-y-y(e*) = 1. Now assuming, for example, that at 
1 TeV r^^ = 10, the 7-ray luminosity of the source below e* = 10 GeV (for 
s = 0.5) should exceed the luminosity of the "TeV photosphere" by an order 
of magnitude. This obviously contradicts the 7-ray observations of Markarian 
421, and especially Markarian 501. Therefore we may conclude that the optical 
depth at 1 TeV, njev = r^-yil TeV) < 10. 

From Eqs.(ll) and (12) we obtain a simple relation between the photo-meson 
cooling time of protons, tp-y, and the optical depth Tuev' 

tp, ~ 1.8 X 10^ Ai3h Sli'' TfT^v ^iV^' s. (13) 

For a steeper spectrum of low-frequency radiation, e.g. a power-law with s — 1, 
and for the same normalization to ruev, the photo- meson cooling time is 
shorter: 

tp^-lO^ AishrfT^v^rg's. (14) 

Even so, remains significantly larger than f*. especially if we take into 
account that for s = 1 the optical depth depends more strongly on energy, 
oc e, and therefore the detection of TeV radiation from Markarian 501 
beyond 10 TeV is an indication that Tuev should be significantly less than 1. 

The source transparency condition for multi-TeV 7-rays implies that in the 
TeV blazars the photo-meson processes proceed on significantly larger time- 
scales compared with t* and tsy (see Fig. 4). In particular, for the power- law 
spectrum of low-frequency radiation with a (most likely) slope s — 0.5 

^d.TxlO^C-^r-v^^rr, (15) 

and 

^ :^ 4 X lO^Aiah B^^, r^^^ . (16) 

Thus, we conclude that for any reasonable assumption concerning the geome- 
try of 7-ray production region, as well as the spectral shape of low-frequency 
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radiation in the blob, the detection of TeV 7-rays from any blazar would im- 
ply low efficiency of the photo-meson processes in the jet, unless the energy 
of protons does not significantly exceed 10^^ eV. In compact 7-ray production 
region(s) of the jet with a typical size it! < 3 x 10^^ At^hdio cm, the proton 
acceleration to such high energies is possible only in the presence of magnetic 
field i? ^ 10 G (see below). At such conditions, however, the synchrotron radi- 
ation becomes a more effective channel for conversion of the kinetic energy of 
accelerated protons to very high energy 7-rays. In principle, the difficulty with 
synchrotron losses could be overcome by adopting a weak {B < IG) magnetic 
field in the blob, but assuming that the EHE protons are accelerated outside 
of the blob, e.g. near the central compact object, and then transported along 
with the jet (Kazanas & Mastichiadis 1999). However, this assumption docs 
not solve the second problem connected with the low efficiency of the photo- 
meson processes in the jet imposed by the transparency condition for TeV 
7-rays. 

A regrettable consequence of this conclusion is the suppressed TeV neutrino 
flux. However, it should be emphasized that this statement concerns only 
the objects seen in TeV 7-rays. Meanwhile, the pair cascades initiated by 
secondary electrons and 7-rays from p7 interactions may still remain a viable 
possibility for other AGN, in particular for the powerful GeV blazars detected 
by EGRET (Mukherjee et al. 1997), as well as for radio-quiet AGN (e.g. Sikora 
et al. 1987), where the radiation density is much higher than in the BL Lac 
objects, and, more importantly, the photo-meson cooling time of EHE protons 
is not constrained by the severe TeV 7-ray transparency condition. 



4 Self regulated synchrotron cutoff 

The energy spectrum of synchrotron radiation depends on the spectrum of 
accelerated protons and the jet's Doppler factor. The high energy cutoff in 
the spectrum of protons is determined by the balance between the particle 
acceleration and cooling times. It is convenient to present the acceleration 
time of particles tacc in the following general form 

t.oo = V(E) rjc = 1.36 X lO" B'^, ri(E) s , (17) 

where Vg = E/(eBj_). The so-called gyro- factor t]{E) > 1 characterizes the 
energy-dependent rate of acceleration. For almost all proposed models of par- 
ticle acceleration in different astrophysical environments, 77 remains a rather 
uncertain model parameter. This is true especially for the small-scale jets of 
blazars, where the nature of the acceleration mechanism itself remains highly 
unknown. On the other hand, any postulation of acceleration of EHE protons 
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in compact 7-ray production regions of the small-scale jets actually implies 
that the parameter 77 should be close to 1, which corresponds (independent of 
a specific mechanisms of acceleration) to the maximum (theoretically possible) 
acceleration rate based on a simple geometrical consideration (Hillas 1984). 
In the case of diffusive shock acceleration in the blazar jets the parameter 77 
typically is expected (Henri et al. 1999) to be larger than 10 (see, however, 
Bednarz & Ostrowski 1996). Therefore, perhaps more effective acceleration 
mechanisms should be invoked for production of EHE protons in small-scale 
AGN jets. An interesting possibihty could be the particle acceleration at the 
annihilation of magnetic fields in the fronts of Poynting flux dominated jets 
(Blandford 1976; Lovelace 1976). It has been argued that this mechanism could 
provide effective acceleration of EHE protons with 77 ~ 1 (Haswell et al. 1992). 

The discussion of particle acceleration mechanisms is beyond the framework of 
this paper, and therefore here the assumption of 77 < 10 should be treated just 
as an working hypothesis. If true, during the characteristic time t* given by 
Eq. (9), the protons could be accelerated in the jet with B ~ 100 G and 5j 10 
up to energies E ~ 10^° cV. At such conditions the losses of highest energy 
protons are dominated by the proton-synchrotron radiation, and therefore the 
cutoff energy Eq is determined by the condition tsy — tacc '■ 



Eo = (3/2)3/4 / 1 ^2^4 ^ 18 ^ 1019 ^-1/2 _ (18) 



Note that the relevant cutoff in the electron spectrum appears much earlier, 
E,,o = {mJm^fEo ~ 5.3 x lO^^ Srof 7?-V2eV. 

Substituting Eq. (18) into Eq. (2) we find that the position of the cutoff 
in the spectrum of the proton-synchrotron radiation is determined by two 
fundamental physical constants, the proton mass nipC^ — 938 MeV, and fine- 
structure constant ojf = 1/137. It depends only on the parameter 77, but not 
on the magnetic field: 

9 

eo = - af^ rripC^ 7]'^ ~ 0.3 7]'^ TeV . (19) 



For the electron-synchrotron radiation, the universal cutoff appears at Sq — 
{mjmp) eJP^ = 9/4 af ^ m^c^ rj-^ ~ 0.160 77-^ GeV. 

Gomparing now the minimum (energy-independent) escape time of protons 
tesc ^ R/c 3.3 X 10^ i?i5^ s, with the synchrotron cooling time of the pro- 
tons responsible for production of 7-rays in the cutoff region, tgy ~ 2.4 x 
10^ BiqJ^t]^/'^ s we find the following condition for formation of the self- 
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regulated synchrotron cutoff : 

5ioo>0.8i?itV/', 



(20) 



where R15 — cm is the source radius in units 10^^ cm. 

If the synchrotron 7-rays are produced in the relativistically moving source, 
the spectral cutoff is shifted towards the TeV domain: 

eo ~ 3 Sio TeV . (21) 

Thus, in the regime of acceleration with 77 < 10, the proton-synchrotron ra- 
diation emitted by a highly magnetized {B ~ 100 G) compact blob with a 
typical size Ri^ ~ 3.2At3h(5io cm and Doppler factor 6j > 10, results in effec- 
tive production of TeV 7-rays. Remarkably, despite possible changes of some 
principal model parameters, first of all the size R and the magnetic field B of 
the production region (e.g. caused by expansion or compression of the blob), 
we should expect a rather stable position of the synchrotron cutoff, provided 
that the parameter rj and the Doppler factor 6^ remain unchanged. Meanwhile 
any change of B and/or R should result in strong (oc B'^R^) variations of the 
absolute flux of synchrotron radiation. This intrinsic feature of the proton- 
synchrotron radiation could explain in a natural way the effect of weak cor- 
relation between the spectral shape and the absolute flux of Markarian 501 
revealed during the extraordinary outburst of the source in 1997 (Aharonian 
et al. 1999a,b). 

Dramatic changes of conditions in the 7-ray production region should lead, in 
fact, to the transition of the source from the regime dominated by synchrotron 
losses to the regime dominated by the particle escape or by adiabatic losses. For 
example, in the case of relativistically expanding blob in which the reduction 
of the magnetic field is faster than B oc at some stage of the source 

evolution the condition represented by Eq.(20) could be violated, and thus 
the spectral shape of th proton-synchrotron radiation would become time- 
dependent. Indeed, in this regime the cutoff energy in the spectrum of protons 
is determined by the balance between tacc and t^sc, which gives 

eo ~ 5.2 BloRl^7j-^5,o TeV . (22) 

Note that Eq. (20) keeps the synchrotron cutoff energy, which is formed in 

the regime dominated by particle escape, always less than the cutoff energy 
formed in the synchrotron-loss-dominated regime. For example, at stages of 
source evolution when i? < 3 x 10^^ cm and magnetic field i? < 10 G, the 
7-ray spectrum is formed in the regime dominated by particle escape, thus 
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e < O.Obr]"^ 5ioTeV, and correspondingly at energies beyond 1 TeV the 7- 
ray spectrum becomes very steep. In the escape-loss- dominated regime the 
cutoff energy, and thus the shape of the entire synchrotron spectrum, strongly 
depend on the magnetic field and the size of the source. In these stages of 
source evolution a significant spectral variability is expected Depending on 
the relationship between the three appropriate time-scales, t^cc, 4y, and t*. wc 
may predict, similar to the electron synchrotron radiation in blazars (Kirk et 
al. 1998), quite different variability patterns such as "soft lags" (when tgy ^ 
t* > tacc), "hard lags" (when 4y ~ ^acc ~ t*), etc. 

Finally, it should be noted that the standard self-regulated synchrotron cutoff 
determined by Eq. (19) corresponds to the scenario when the proton accelera- 
tion and 7-ray production take place in the same localized region. In the case 
of external injection of accelerated protons into the magnetized cloud, the en- 
ergy of synchrotron radiation is determined by the magnetic field in the 7-ray 
production region and the maximum energy of injected protons, therefore it 
could exceed the 'self-regulated synchrotron cutoff. 



5 Implications to Markarian 501 and Markarian 421 

The HEGRA and Whipple observations of Markarian 501 during the extraor- 
dinary outburst in 1997 revealed that despite spectacular (up to factor of 
10 or more) flux variations, the shapes of the daily 7-ray spectra remained 
unchanged^ throughout the entire state of high activity (Aharonian et al. 
1999a,b; Quinn et al. 1999). To some extent this could be true also for Markar- 
ian 421 (Aharonian et al. 1999c; Krcnnrich et al. 1999), although the uncer- 
tainties in the daily TeV spectra are too large for a firm conclusion. 

Markarian 421 and Markarian 501 exhibit substantially different 7-ray spectra. 
The time-averaged spectra of Markarian 501 in 1997 (Aharonian et al. 1999b) 

and Markarian 421 (Aharonian et al. 1999c) as measured by the HEGRA 
system of telescopes are shown in Fig. 5. The differential flux of Markarian 
501 from 500 GeV to 24 TeV is described by a power- law with an exponential 
cutoff: dN/de oc exp (— e/eo) with F ~ 2 and eo ^ 6.2 TeV (Aharonian et 
al., 1999b). Remarkably, a quite similar spectrum was detected by the HEGRA 
collaboration during a short, but strong outburst of Markarian 501 in June 
1998 (Sambruna et al. 2000). In the TeV regime, Markarian 421 has a steeper 
spectrum. In the interval from 500 GeV to 7 TeV it could be approximated by 

^ In fact, the CAT group has found a tendency of the spectrum of Markarian 501 
to become somewhat harder during the flares (Djannati-Atai et al. 1999). However, 
the effect is not very strong and could not change the general conclusion about 
essentially stable spectral shape of TeV radiation in the high state. 



15 




Photon Energy, eV 



Fig. 5. The proton-synchrotron radiation fits to the TeV spectra of Markarian 501 
and Markarian 421. The TeV data of Markarian 501 correspond to the high-state 
spectrum of the source as measured by HEGRA in 1997 (Aharonian ct al. 1999b). 
The TeV data of Markarian 421 correspond to the average HEGRA spectrum based 
on the measurements in 1997-1998 (Aharonian et al. 1999c). The dashed and sohd 
curves correspond to the spectra of proton-synchrotron radiation before and after 
correction for intergalactic extinction, respectively. A "low-CBR" model is assumed 
(solid curve in Fig. 6). The data of Markarian 501 are fitted assuming Op = 2 
and eo = 1.3 TeV, and the data of Markarian 421 are fitted assuming ap = 2 and 
eo = 0.26 TeV. 

a pure power-law with a photon index F ~ 3 (Aharonian et al. 1999c; Piron 
et al. 2000), or by a power-law with F = 2.5 ± 0.4 and exponential cutoff at 
eo = 2.8lo:9TeV (Aharonian et al. 1999c). 

5. 1 Extragalactic extinction of gamma rays 

The spectra of TeV radiation observed from distant {d > 100 Mpc) extragalac- 
tic objects suffer essential deformation during the passage through the inter- 
galactic medium caused by energy-dependent absorption of primary 7-rays 
at interactions with the diffuse extragalactic background radiation (Gould & 
Schreder 1966; Stecker et al. 1992, Vassiliev 2000). For Markarian 421 and 
Markarian 501 at distances 150 Mpc and 170 Mpc, respectively, this effect 
could be significant already at sub- TeV energies (Guy et al. 2000), and be- 
comes especially strong at energies above 10 TeV, where the optical depth 
r-y-y most hkely significantly exceeds 1 (Goppi & Aharonian 1999b). The lack 
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of a relevant broad-band information about the cosmic background radiation 
(CBR) from sub-micron to sub-mm wavelengths introduces an ambiguity in 
the interpretation of the observed 7-ray spectra. Nevertheless, the recent re- 
ports about detection of CBR at near infrared (2.2 /xm and 3.5 //m - Dwek & 
Arendt 1998; Gorjian et al. 2000), mid infrared (15 //m - Elbaz et al. 1999), 
and far infrared (140 //m and 240 /xm - Schlegel ct al. 1998; Hauser et al. 
1998) bands, allow significant reduction of the uncertainty in the intergalactic 
extinction. Yet, these measurements alone are not sufficient for quantitative 
study of the intergalactic extinction of 7-rays, n{E) = exp [— r^^(e)], in the 
observed energy range from 0.5 to 20 TeV. Apparently, a suitable theoretical 
model of CBR is needed which should reproduce properly the shape of both 
"stellar" and "dust" components of the background radiation. Note that most 
of cosmological models give rather similar shapes of CBR with two distinct 
bumps at 1-2 /im and 100 -200 /im and a valley at 10-20 /im (see, e.g., Dwek 
et al. 1998; Primack et al. 1999). 




A (/um) 

Fig. 6. The Cosmic Background Radiation (CBR). The solid curve corresponds 

to the "low-CBR", and the dashed curve corresponds to the "high-CBR" model 
predictions (see the text). The dotted curve shows the flux of the 2.7 K microwave 
background radiation. The tentative detections of CBR at 2.2/xm and 3.5/Ltm (filled 
circles) are from Gorjian et al. (2000), and at 15/im (the star) from Elbaz et al. 
(1999); the fluxes at 140^m and 240;um (filled triangles) are from Hauser et al. 
(1998). The other upper/lower CBR fiux limits (open symbols) are taken from the 
recent compilation of CBR by Dwek et al. (1998). 

Below I will use the so-called LCDM model of Primack et al. (1999), but 
the absolute fluxes of both "stellar" and "dust" components of radiation will 
be allowed to vary within a factor of 2 or so. The prediction for the CBR 
flux by the LCDM model (hereafter "low-CBR") and re-scaled LCDM model. 
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assuming 50 per cent higher flux of the "stellar" component, and twice higher 
flux of the "dust" component (hereafter "high-CBR"), are presented in Fig. 6. 
For comparison, the recently reported CBR flux measurements or upper/lower 
limits from UV to FIR are also shown. Despite large statistical and systematic 
uncertainties, the actual spectrum of CBR most probably does not deviate 
significantly from the two simplified model predictions shown in Fig. 6. 



5.2 Fitting the TeV spectra of Markarian 501 and Markarian 421 

The spectrum of the proton- synchrotron radiation corrected for intergalactic 
extinction can satisfactorily fit, for a reasonable combinations of a limited 
number of model parameters, the observed TeV fluxes of both Markarian 501 
and Markarian 421. This is demonstrated in Fig. 5, where the average TeV 
fluxes of both objects are shown together with the model predictions. The 
theoretical spectra, after corrections for the intergalactic extinction assuming 
the "low CBR" , are normalized to the measured fluxes at 1 TeV. This deter- 
mines, for the given magnetic held, the required total energy in accelerated 
protons. Assuming that the accelerated protons have an energy distribution 
represented by Eq. (7), the observed spectra of TeV emission can be fitted 
by the proton-synchrotron radiation for the following combinations of model 
parameters: Op = 2 and eo = 1.3 TeV for Markarian 501, and ftp = 2 and 
eo = 0.26 TeV for Markarian 421. Assuming now that the 7-rays are pro- 
duced in the synchrotron-loss-dominated regime, from Eq. (21) one can easily 
estimate the ratio p — 5io/r) which is the most relevant parameter for deter- 
mination of the position of the self-regulated synchrotron cutoff - p — 0.43 
and 0.09 for Markarian 501 and Markarian 421, respectively. 

The predicted for Markarian 501 proton-synchrotron radiation fiux passes be- 
low the HEGRA point at 24 TeV (see Fig. 5). This discrepancy could be 
removed, in principle, by assuming a different (steeper) shape of the spectrum 
of the "dust" component compared with predictions of current CBR models 
(Coppi & Aharonian 1999b), and/or assuming a flatter spectrum of protons 
in the region of the cutoff Eq compared to the spectrum given by Eq. (7). 
However, because of the inadequate statistical significance of the 7-ray data 
above 16 TeV (Aharonian et al. 1999b), the last point on Fig. 5 should be 
considered a fiux upper limit rather than a positive detection. 

If one treats the TeV observations independently of data obtained in other en- 
ergy domains, the parameters chosen to fit the TeV spectra cannot be unique 

in a sense that some other combinations of Op and cq could still provide rea- 
sonable spectral fits. This is demonstrated in Fig. 7 where the measured fluxes 
of Markarian 501 are shown together with the predicted spectra of the proton- 
synchrotron radiation assuming "high-CBR" model for the background radi- 
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Fig. 7. The SED of Markarian 501. The thin and heavy curves correspond to the 
spectra of proton-synchrotron radiation before and after correction for the inter- 
galactic extinction, respectively. The "high-CBR" model is assumed (dashed curve 
on Fig. 6). The data of Markarian 501 are fitted by parameters Op = 2 and 



eo = 1.8 TeV (sohd curve) 



Or 



1.5 and eg = 1.2 TeV (dashed curve), and = 2.4 



and eo = 3.5 TeV (dot-dashed curve). 



ation. It is seen that the detected fluxes could be equally well fitted by syn- 
chrotron radiation for 3 different combinations: (i) ttp = 2 and eo = 1.8 TeV 
(or p ~ 0.6, if the energy losses of protons are dominated by the synchrotron 
radiation), (ii) = 1.5 and eo = 1.2 TeV (p ~ 0.4), and ap = 2.4 and 
Co = 3.6 TeV (p ~ 1.2). The impact of the uncertainty in the CBR flux can 
be estimated by comparing the solid curves in Fig. 5 and Fig. 7. While in 
both cases the power-law index of the proton spectrum is assumed ftp = 2, 
the "high-CBR" requires somewhat (~ 40 per cent) higher synchrotron cutoff 
compared with the case of "low-CBR". Also, for the "high-CBR" the dis- 
crepancy of the theoretical and measured fluxes above 16 TeV becomes more 
apparent. 

An important constraint on the possible "Op/eo" combination could be pro- 
vided by 7-ray data obtained at low, MeV/GeV energies. In particular, the 
combination of parameters used for the fit of the spectrum of Markarian 421 
in Fig. 5, predicts an energy fiux at 100 MeV of about 2 x 10"^-*^ erg/s, i.e. 
factor of 3 less than the EGRET fiux averaged over the 1992-1996 period of 
observations (Sreekumar et al. 1996). Therefore, assuming that the observed 
MeV/GeV fiuxes are associated with the proton-synchrotron radiation as well, 
the best fit of the data from 100 MeV to 10 TeV could be obtained adopting 
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a steeper proton spectrum with ccp = 2.4 (see Fig. 8). Although the data ob- 
tained at TeV and MeV/GeV energies correspond to different epochs, the lack 
of evidence of strong spectral variability in both energy domains (Aharonian 
et al. 1999c; Sreekumar at al. 1996) makes such a fit quite meaningful. 
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Fig. 8. The proton-synchrotron radiation of Markarian 501 and Markarian 421 from 
MeV/GeV to TeV energies. All spectra are normalized to the observed fluxes at 
1 TeV after correction for the intergalactic extinction assuming the "low-CBR" 
model. The data of Markarian 421 are fitted by ap = 2.4 and eo = 1 TeV. The data 
of Markarian 501 are fitted by three different combinations of ap and cq: Op = 2 
and eo = 1.3 TeV (dashed curve); ap = 2.6 and eg = 4TcV (dot-dashed curve); 
ap = 1.5 and eg = 0.75 TeV (dotted curve). The TeV data of Markarian 501 and 
Markarian 421 are the same as in Fig 2. The filled square at 100 MeV corresponds 
to the flux upper limit set by EGRET during the period April 9-15 (1997) when 
Markarian 501 was in a very high state (Catanese et al. 1997). The zone of the 
low-energy 7-ray fluxes of Markarian 421 shown from 50 MeV to 5 GeV corresponds 
to the spectrum of the source averaged over the Phase-1 period of the EGERET 
observations (Sreekumar et al. 1996). 

The observations of Markarian 501 by EGRET during the extraordinary active 
state in April 1997 resulted only in the flux upper hmit at 100 MeV (Catanese 
et al. 1997). This upper limit, combined with the TeV data, constrains the 
power-law index of the proton spectrum, ap < 2.6. For «p = 2.6 the position 
of the synchrotron cutoff should appear at high energies, cq = 4 TeV, in order 
to match the TeV spectrum (see Fig. 8). For the very flat proton spectrum 
with CKp = 1.5, the synchrotron cutoff occurs below 1 TeV, eo ^ 0.75 TeV, 
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otherwise the calculated spectrum becomes harder than the observed TeV 
spectrum. And finally, for the "standard" proton spectrum with = 2, the 
best fit value for the synchrotron cutoff is provided by eg ~ 1.3 TeV, or p ~ 
0.43 if the TeV spectrum of Markarian 501 is formed in the synchrotron- 
loss-dominated regime. Because the Doppler factor of the small-scale jet in 
Markarian 501 is believed to be close to 10, a conclusion could be drawn that 
during the entire 1997 outburst the acceleration of particles in Markarian 501 
took place in the regime of maximum acceleration rate [t] ~ 2). If so, this 
explains in a rather natural way the essentially time-independent shape of 
the TeV spectrum observed during the strong flares in 1997 (Aharonian et al. 
1999a,b). 

Meanwhile, a significant drop in the acceleration rate (and/or the strength of 
the magnetic field) should lead to dramatic changes of the spectrum of 7-rays. 
This effect perhaps can explain the recent report of detection by EGRET of 
a very fiat 7-ray fiux at energies from 50 MeV to 5 GeV during the period 
of relatively low state of Markarian 501 in March 1996. A rather surprising 
aspect of this report is that the energy fiux of 7-rays at several GeV was sig- 
nificantly larger than the X-ray and TeV fiuxes measured approximately at 
the same time by ASCA and Whipple instruments (Kataoka et al. 1999). In 
the framework of the "proton-synchrotron radiation" hypothesis, the EGRET 
and Whipple fiuxes could be explained with model parameters = 2 and 
eo = 0.02 TeV, or, if the proton acceleration takes place in the regime domi- 
nated by synchrotron losses, p ~ 0.007 (Fig. 9). The latter implies dramatic 
reduction of the acceleration efficiency {rj ~ 1405io) compared with the high 
state, provided that the jet's Doppler factor remains more or less constant. 
It is necessary to notice that this interpretation which adopts very small p 
parameter, but tacitly assumes that the proton acceleration takes place in the 
synchrotron-loss-dominated regime, requires, as it follows from Eq. (20), an 
extremely large magnetic field, B > AOOR^^ S^q G. A more realistic scenario 
seems the transition of the source from the synchrotron-loss- dominated regime 
to the particle escape-loss- dominated regime, for example, due to a possible 
reduction of the magnetic field caused by the expansion of the blob. The sig- 
nificant drop of the magnetic field, e.g. from B ~ 100 G to i? ~ 10 G would 
result in an earlier synchrotron cut-off. 

The signature of such a low state is the shift of the position of the peak of 
SED down to 100 GeV or below, with a hard spectrum of GeV 7-rays, and 
very steep spectrum at TeV energies. While such a hard spectrum of GeV 
radiation is indeed observed, the small TeV photon statistics do not provide 
an adequate information about the spectrum at highest energies. Contrary to 
the synchrotron-loss- dominated regime, the spectrum of 7-rays formed in the 
escape-loss-dominated regime depends on parameters characterizing the pro- 
duction region, in particular the source size and the magnetic field. Therefore, 
in this stage of evolution, the source should show significant variations of the 
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Fig. 9. The proton-synchrotron radiation of Markarian 501 in a high and low states. 
The dashed and soHd curves correspond to the spectra of proton-synchrotron ra- 
diation before and after correction for intergalactic extinction, respectively. The 
"low-CBR" model is assumed. The data of Markarian 501 in the high state are the 
same as in Fig. 2. The broad-band 7-ray data in the low state are obtained by the 
Whipple (open triangle) and EGRET (filled triangles) groups during the multiwave- 
length campaign in March 1996 (Kataoka et al. 1999). An "archival" upper limit 
on the 100 MeV flux based on the long-term observations of the source during the 
phase I period of EGRET is also shown (star) . The high state of the source is fitted 
by Q!p = 2 and eo = 1.3 TeV. The low state of the source is fitted by ftp = 2 and 
eo = 0.02TeV. 

spectral shape of TeV 7-ray emission. 

The TeV spectrum of Markarian 421 is significantly steeper than the spectrum 
of Markarian 501. Another distinct feature of Markarian 421 is the TeV fiux 
variability on extremely short, At ~ 15 minute timescales (Gaidos et al. 1996). 
Assuming that the radiation is formed in the synchrotron-loss-dominated 
regime, the steep 7-ray spectrum of Markarian 421 with eo ~ 0.26 TeV (see 
Fig. 5) requires relatively slow acceleration rate, ~ 12(5io. On the other hand, 
in this regime the radiative cooling time of protons responsible for production 
of 7-rays in the cutoff region, tgy — 2.4 x 10^ B-^qJ'^i]^^'^ s, could match the 
observed variability (in the frame of the jet) 9 x lO^^io Sj if -Bioo ^ '^V^^'^^io^^ ■ 
Therefore, i?ioo > 4.5 S^q^^. Thus, if the TeV radiation of Markarian 421 is 
formed in the synchrotron-loss-dominated regime, the magnetic field should 
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significantly exceed 100 G, which seems a rather extreme assumption. 

In a relatively moderate (and perhaps more realistic) magnetic field of about 
100 G or less, both the observed rapid variability and the steep TeV spectrum 
of Markarian 421 could be interpreted, assuming that the radiation is formed 
in the escape-loss-dominated regime. Indeed, in the compact 7-ray production 
region with i? < c At 5j ~ 2.7 x lO^^^io cm, the position of the synchrotron 
cutoff is estimated eo — 0.38 i^j^go'^io^^^ '^^^ (^^^ ^2). For a reasonable 
combination of parameters Sioo ~ 1, (^10 ~ 1 and ~ 1, this estimate agrees 
quite well with the allowed range of the synchrotron cutoff between 0.25 and 1 
TeV as derived from the TeV spectrum of Markarian 421 (Fig. 5 and Fig. 8). 

5.3 Synchrotron radiation from the proton distribution with sharp pile-up 

The spectra of the proton-synchrotron radiation discussed above are calculated 
under the assumption of a "standard " (power-law with exponential cutoff) 
proton spectrum represented by Eq. (7). Actually, the spectra of accelerated 
particles at some circumstances could have more complicated and exotic forms, 
e.g. they could contain pile-ups preceding the spectral cutoffs. An example of 
a SED of synchrotron radiation corresponding to the proton spectrum with 
sharp pile- up and abrupt cutoff is shown in Fig. 10. The shape of this spectrum 
(dashed curve) essentially differs from the measured spectrum of Markarian 
501. Moreover, for any reasonable spectral distribution of the GBR, the de- 
formation of the initial 7-ray spectrum caused by the intergalactic extinction, 
hardly could bring it closer to the observed one. The desirable shape, how- 
ever, could be achieved, if we assume that the deformation of the primary 7-ray 
spectrum takes place in two steps, e.g. presuming that the primary radiation 
suffers significant absorption twice - firstly in the vicinity of the source, and 
afterwards - in the intergalactic medium. The synchrotron radiation of the jet 
itself seems to be a natural photon field to serve as an internal absorber for the 
TeV radiation. Actually, such absorption is even unavoidable if the Doppler 
factor of the jet Sj < 10. However, the pronounced peak in the production spec- 
trum of radiation shown in Fig. 10 (dashed curve), requires rather selective 
absorption with strongest impact in the region of the 7-ray peak, which can- 
not be provided by the broad and flat spectrum of low-frequency synchrotron 
radiation of the jet. Because 77 e+e^ cross-section peaks near the thresh- 
old, there is approximately a one-to-one correspondence between the energies 
of the absorbed 7-ray e* and the field photon hu: e* ~ 1 (/iz//l eV)~^ TeV. 
Therefore, for a selective absorption of 7-rays from the region of the maxi- 
mum of the production spectrum, the target photon field should have a quite 
narrow distribution, Au/u ~ 1. The dot-dashed curve in Fig. 10 demonstrates 
the effect of such absorption in a monochromatic photon field with energy 
hu — 0.15 eV, and assuming that r-y-y = 1 at e ~ 7 TeV. A more realistic. 
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Fig. 10. The proton-synchrotron radiation of Markarian 501 calculated for the pro- 
ton spectrum with a sharp pile-up (curve 3 in Fig. 1). Dashed curve - the production 
spectrum; dot-dashed curve - 7-ray spectrum after internal absorption; solid curve 
- 7-ray spectrum after internal and intergalactic absorptions. The following param- 
eters arc assumed: cq = 7.8 TeV (or p = 2.6), the mean energy of target photons 
responsible for internal absorption hu = 0.15 eV, optical depth r^^ = 1 at e = 7 TeV. 
For the extragalactic extinction of 7-rays the "high-CBR" model is assumed. 



e.g. thermal (Plankian) distribution of the target photons with temperature 
kT ~ 1/3 /iz/ ~ 0.05 eV gives a similar result, with a bit smoother spectrum of 
absorbed 7-rays at 2 TeV. The optical depth r^^ ~ 1 caused by mid infrared 
(A ~ 10/xm) photons, could be provided by different kind of sources in the 
vicinity of the central black hole (Celotti et al. 1998). 



The TeV 7-rays, after suffering additional absorption in the intergalactic pho- 
ton fields, will eventually arrive with "double-deformed" spectrum shown in 
Fig. 10 (sohd curve). The calculated synchrotron spectrum, after corrections 
for internal and intergalactic absorption and normalized to the observed flux 
at 1 TeV, quite satisfactorily fits the SED of Markarian 501 in the high state, 
taking into account significant (up to factor of 2) systematic uncertainties of 
the HEGRA fiux measurements at energies below 1 TeV and above 15 TeV 
(Aharonian et al., 1999b). 
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6 Discussion 



Although there is httle doubt that the bulk of the highly variable X-ray emis- 
sion of the BL Lac objects has synchrotron origin and is produced in the 
relativistic jets pointed to the observer (Urry & Padovani 1995), the TeV radi- 
ation of Markarian 421 and Markarian 501 is the only model- independent and 
unambiguous indicator of acceleration of ultra-relativistic particles in these 
small-scale (sub-pc) jets. In the currently popular synchrotron-self-Compton 
(SSC) models of BL Lac objects both the X-ray and 7-ray components are at- 
tributed to the radiation of directly accelerated TeV electrons. For a certain, 
physically well justified combination of parameters (the blob size, magnetic 
field, Doppler factor, etc.), these models (sec, e.g., Inoue & Takahara 1996, 
Mastichiadis & Kirk 1997, Pian et al. 1998, Bednarek & Protheroe 1999, Coppi 
& Aharonian 1999a, Wagner et al. 1999; Kataoka et al. 1999, Maraschi et al. 
1999, Takahashi et al. 1999, Krawczynski et al. 2000) give quite satisfactory ex- 
planation of the observed spectral and temporal characteristics of non-thermal 
radiation over more than ten decades of frequencies from 10^^ Hz to 10^'' Hz. 
In the SSC models the strength of the magnetic field in the 7-ray emitting 
regions typically is less than 1 G, the most likely value being close to 0.1 G, 
which directly follows from the comparable energy fluxes released in X-rays 
and 7-rays during the TeV flares of Markarian 421 and Markarian 501 (e.g. 
Tavecchio et al. 1998). At such a low magnetic field the maximum energy 
of accelerated protons cannot exceed -Ep^max ~ lO^^Atsh^io {B/lG)ri^^cV, 
and therefore the contribution of protons in 7-ray production through both 
the photo-meson and the synchrotron channels is negligible; the production of 
TeV 7-rays in the jet is strongly dominated by the inverse Compton scattering 
of electrons directly accelerated up to energies ~ 10 TeV. 

The increase of the magnetic field leads to reduction (oc B~^) of the fiux of 
the inverse Compton 7-rays; for i? > 1 G the bulk of energy of accelerated 
electrons is channeled into the synchrotron radiation. Also, the increase of the 
magnetic field shifts (oc B) the synchrotron peak to higher energies. But for 
any reasonable Doppler factor of the jet, the observed TeV fluxes of Markarian 
421 and Markarian 501 cannot be explained by the electron synchrotron radi- 
ation because of the self-regulated synchrotron cutoff at Cmax — 1-6 5io GeV, 
which inevitably appears if the electron acceleration and 7-ray production take 
place in the same region of the jet. Formally, we may avoid this limit assuming 
that the regions of the electron acceleration and the 7-ray production are sep- 
arated, i.e. the electrons are accelerated up to energies ~ 10^ TeV in a region 
with rather small magnetic field, B < 0.01 G, but release all their energy in a 
form of GeV/TeV synchrotron radiation later, after entering the region(s) of 
strongly compressed/amplified magnetic field, i? > 10 G. This scenario seems, 
however, rather artificial. 
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Large magnetic fields, B ~ 100 G or so, coupled with effective acceleration of 
protons at the maximum rate, may create very favorable conditions for TeV 7- 
ray production. Indeed, at such conditions the synchrotron cooling of protons 
not only well dominates over other radiative and non-radiative losses (see 
Fig. 4), but also provides good fits to the observed spectra of TeV radiation 
of Markarian 421 and Markarian 501. Remarkably, if the proton acceleration 
takes place in the synchrotron-loss- dominated regime, the spectral shape of 
radiation depends only on the power-law index of accelerated protons, ap and 
the parameter p — Siq/t), but not on the magnetic field and the size of the 7-ray 
emitting blobs, which generally endure significant time-evolution. Meanwhile, 
any change in the size and/or magnetic field of the evolving blobs should 
lead to significant variation in the absolute flux of 7-rays. This effect could 
give a natural explanation for one of the remarkable features of strong flares 
of Markarian 501 (and, perhaps, also Markarian 421) - the essentially stable 
spectral shape of TeV radiation despite strong variation of the absolute flux 
observed on time scales less than 1 day. During the strong flares of Markarian 
501 with a flat spectrum around 1 TeV, the parameter p = 610 /r] ~ 1. The 
exact value of p depends on the level of the diffuse extragalactic background 
and on the power-law index of accelerated protons. The uncertainty in this 
parameter, despite signiflcant uncertainties in the intergalactic extinction of 
TeV radiation, as well as the lack of adequate information about the 7-ray 
fluxes at MeV/GeV energies, does not exceed a factor of three or so. Since 
the Doppler factor in Markarian 501 is believed to be within 10 to 30, the 
acceleration rate should be very high, rj = p^^ (5io ~ 1 — 3. The unusually high 
GeV flux observed from Markarian 501 during the multiwavelength campaign 
in March 1996, when the source was in a very low TeV state, could be explained 
by a transition of the source from the synchrotron-loss- dominated regime to the 
escape-loss- dominated regime, caused, for example, by dramatic drop of the 
magnetic field. In this regime we should expect steeper TeV 7-ray spectra with 
a slope depending on the magnetic field and the size of the source. Therefore, 
an important test of this hypothesis could be detection of spectral variability 
of 7-radiation in a low state of the source. 

The significantly steeper TeV spectrum of Markarian 421 requires early syn- 
chrotron cutoff, eo ~ 0.25 — 1 TeV. This could be interpreted as a result of 
formation of the synchrotron spectrum in the regime dominated by the es- 
cape of protons from the acceleration region. For a magnetic field of about 
100 G this would require a compact source with R < 10^^ cm. Interestingly, 
such a small linear size of the 7-ray production region is supported indepen- 
dently by the observed dramatic variations of the TeV fiux of Markarian 421 
on timescales ~ 15 minutes. 

Within the model of proton synchrotron radiation, we may expect also syn- 
chrotron radiation produced by directly accelerated electrons - the counterparts 
of EHE protons. If the particle acceleration takes place in the synchrotron- 
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loss- dominated regime, the self-regulated synchrotron cutoff of this component 
depends only on the parameter p = Siq/i], namely eg — 1.6pGeV, thus it cor- 
relates with the position of the cutoff in the proton synchrotron spectrum at 
eo ^ 3p TeV. The ratio of the energy fluxes of these two components is deter- 
mined simply by the ratio of non-thermal energy channeled into the acceler- 
ated protons and electrons, Wp/We- In a strong magnetic field of about 100 G, 
the synchrotron cooling time of electrons is shorter, almost at all rclativistic 
energies, than the typical dynamical (e.g. light- crossing) times. This results in 
a well-established steady-state spectrum of electrons dN/dE oc £'~(°o+i), pro- 
vided that the power- law index of acceleration spectrum ckq > 1. Consequently, 
a pure power law synchrotron spectrum with a photon index (ao + 2)/2 would 
be formed. In particular, for ao = 2, we should expect a fiat synchrotron SED 
{I'Si, — const) from the optical/UV wavelengths to MeV/GeV 7-rays. The 
broad-band spectra of both Markarian 421 and Markarian 501 do not agree 
with such a pure power-law behavior; in fact, the SED of both objects show 
pronounced synchrotron X-ray peaks. Therefore this (theoretically possible) 
population of directly accelerated electrons - counterparts of the EHE protons 
- cannot be responsible for the bulk of X-ray emission. 

The latter could be referred to electrons, produced, most probably, in a dif- 
ferent way, and/or in other region(s) of the jet. The X-ray light curves of 
both Markarian 421 (e.g. Takahashi et al. 1999) and Markarian 501 (e.g. Sam- 
bruna et al. 2000) show so-called "soft" and "hard" lags. A possible inter- 
pretation of this effect in terms of competing acceleration, radiative cool- 
ing, and escape timescales of synchrotron-emitting electrons (Takahashi et 
al. 1996; Kirk et al. 1998), would require that all these timescales are com- 
parable with the light crossing time, t = R/c ~ 10^ — 10^ s. The cool- 
ing time of an electron responsible for a synchrotron photon of energy e is 
t§ ~ 1.5 X 10-'^(5/lG)-3/2 (e/lkeV)i/2s_ Therefore in the X-ray production 
region the magnetic field cannot, independent of specific model assumptions, 
significantly exceed 0.1 G. Thus, the hypothesis of proton-synchrotron origin of 
TeV radiation implies that the production regions of TeV 7-rays {B ~ 100 G) 
and and synchrotron X-rays (5 ~ 0.1 G) are essentially different. Although in 
any reasonable scenario we may expect a non-negligible correlation between 
the X-ray and TeV 7-ray fluxes, the spatial separation of the X-ray and TeV 
7-ray production regions does not allow definite predictions for such a corre- 
lation. 

In the proton-synchrotron model of TeV radiation of BL Lac objects, two 
more components of X-radiation are expected. Firstly, in the field of about 
100 G the accelerated protons of energy E 10^^ eV themselves produce syn- 
chrotron X-rays. However, the contribution of this component to the observed 
X-ray flux is negligible. A much more prolific channel for X-ray production 
connected (indirectly) with the EHE protons, can be provided by electrons 
of non-acceleration origin, namely by secondary electrons produced at inter- 
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actions of the primary TeV 7-rays with the ambient low-frequency radiation. 
Actually, for a spectrum of EHE protons containing sharp spectral pile-up, we 
must assume an essential (r^^ ~ 1) internal absorption of 7-rays in order to 
match the observed spectrum of Markarian 501 (see Fig. 10). The appearance 
of secondary electrons in the jet results in production of a hard synchrotron 
X-ray component. Apparently, for an optical depth t^^ ~ 1, the luminosity 
of this component would be comparable to to the luminosity of their "grand- 
parents" - TeV 7-rays. The photo-produced electrons have a rather specific, 
significantly different from the directly accelerated particles, shape. For ex- 
ample, the spectrum of electrons produced at interactions of high energy 7- 
rays with a photon index F and field photons with a narrow (e.g. Planckian) 
spectral distribution with a characteristic energy hu, has the following char- 
acteristic form: starting from the minimum (allowed by kinematics) energy at 
-E* = mlc^/Ahu, the electron spectrum sharply rises reaching the maximum 
at ~ 2AE^ ~ 0.15(7w/leV)'^ TeV, and then a,t E :$> it decreases 
as q± oc E~^^'^^^laE. Within an accuracy better than 20 per cent, the spec- 
trum of secondary pairs could be approximated in a simple analytical form 
(Aharonian & Atoyan 1991): 

where x = E/E^ > 1, /(F) = (1.11 - 1.60F + I.ITF^), and J q±{E) dE ^ 2 
(two electrons per interaction). 

The intensive synchrotron losses in the strong magnetic field B ~ 100 G 
quickly establish a steady-state spectrum of electrons proportional to E~^ 
below and approximately as £'^(^+2) (if we ignore the weak logarithmic 
term) above E^. Correspondingly, the synchrotron spectrum of the secondary 
pair-produced electrons is characterized by a smooth transition, through the 
energy around = 120Sioo(/ii^/l eV)-^ keV, from s = 1.5 to s f« (F + 3)/2. 
Assuming, for example, that the TeV 7-rays with a photon index F ~ 1.4 
are absorbed in an external infrared photon field with characteristic energy 
in the jet frame hp = 1.55ioeV (like in Fig. 10), one should expect in the 
observer frame a hard X-ray component of radiation with a SED i/Su oc e°'^ 
and uS^, oc e"^ '^ below and above €5 ~ 500Sioo Siq keV, respectively. Finally, 
at energies above several McV the spectrum becomes very steep due to the 
cutoff in the spectrum of primary (proton-synchrotron) 7-rays. We may spec- 
ulate that such a component of hard X-rays was observed during the strong 
TeV flares of Markarian 501, in particular in April 1997 (Plan et al. 1998). 
Obviously, this component of radiation should strongly correlate with the flux 
of TcV emission. However, because the synchrotron radiation produced by 
directly accelerated electrons significantly contributes to the observed X-ray 
fiux as well, the TeV/X-ray correlation may have a rather complicated and 
non-standard behavior. The absence (or suppression) of such a hard X-ray 
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component in the spectrum of Markarian 421 could be explained, within the 
framework of this model, by the steep spectrum of TeV 7-rays caused by 
an early synchrotron cut-off. The quantitative study of this question within a 
detailed time-dependent treatment of the problem will be discussed elsewhere. 

In BL Lac objects like Markarian 501, the photo- meson processes do not play, 
most probably, dominant role in the production of high energy 7-rays. Indeed, 
the TeV 7-ray transparency condition puts a robust lower limit on the charac- 
teristic time of this process, tp-y ~ 10^ s which is almost 3 orders of magnitude 
larger than the characteristic synchrotron cooling time of protons. The severe 
synchrotron losses of protons cannot be avoided since the requirement of ex- 
tremely high energy protons in the blob - a key assumption in the PIC model 
- automatically implies a large magnetic field of an order of 100 G. On the 
other hand, the characteristic photo-meson cooling time cannot be arbitrary 
reduced, e.g. by assuming extremely high density of radiation in the blob. For 
any reasonable spectrum of the latter, and any reasonable geometry of the 7- 
ray production region, this would lead to an unacceptably large optical depth 
for TeV 7-rays. Therefore, one may conclude that the very fact of observations 
of TeV 7-rays from Markarian 501 and Markarian 421 (in fact, from any TeV 
blazar) almost rules out the "photomeson" origin of the bulk of the observed 
TeV radiation. This also imphes very low TeV neutrino fluxes. 

The hypothesis of proton-synchrotron origin of TeV 7-fiares from Markarian 
421 and Markarian 501 requires a very large amount of energy contained in 
the form of magnetic field 

Wb ~ ^R' = 5.6 X 10^^ ^tlA Bfoo erg ■ (24) 

For comparison, the energy of the magnetic field allowed by the SSC mod- 
els {B ~ 0.1 G) is less by six orders of magnitude. On the other hand, the 
estimates of the kinetic energy in > 10^^ eV protons (within the proton- 
synchrotron radiation model) and > 10^^ eV electrons (within the SSC model) 
are comparable, if we assume that in both models the 7-rays are produced with 
high efficiency. Indeed, in this case the inverse Compton cooling time of elec- 
trons and the synchrotron cooling time of protons responsible for TeV emission 
in the SSC and proton-synchrotron models, respectively, are equal or less than 
the observed variability time-scale t* ~ 1.08 x 10^ Atah^io s, thus 

H/(ssc) _ ^.^PSR) < 4^^2 ^* ^ ^ ^ ^ ^q46 Atsi,S^i erg . (25) 

In this estimate the average flux of Markarian 501 in 1997, corrected for extra- 
galactic extinction for "low-CBR" model is assumed. Higher GBR fluxes, as 
well as the possible internal absorption of 7-rays, would increase this estimate 
by a factor of 3 or so. If the spectrum of accelerated particles extends, e.g. as 
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down to energies E ~ mc^, the above estimate could be increased by an 
order of magnitude. In fact, this simphfied estimate only reflects the average 
level of the energy content of accelerated particles during the high state of 
Markarian 501 in 1997. For the strongest flares with duration < 1 day, the en- 
ergy in accelerated protons could be several times larger than it follows from 
Eq. (25). But in any case, these uncertainties do not prevent us to conclude 
that in the proton-synchrotron model we deal with a highly magnetized con- 
densation of 7-ray emitting clouds of EHE protons, where the magnetic pres- 
sure dominates over the pressure of relativistic protons ^ . The total magnetic 
energy given by Eq. (24) of a single blob formally is sufficient to supply the 
necessary energy in TeV 7-rays released during the entire high activity state of 
Markarian 501 in 1997. Thus speculating that an effective acceleration mech- 
anism operates in such magnetized condensations, stimulated presumably by 
interactions of the latter with the surrounding plasma, a single (or a few) en- 
ergetic blob(s) ejected from the central source towards the observer could, in 
principle, explain the extraordinary high TeV flux of Markarian 501 in 1997. 

In the SSC models the situation is exactly opposite - the pressure of relativistic 
particles is significantly larger than the pressure of the magnetic field. Even so, 
the total energy contained in the relativistic electrons is sufficient to support 
the observed X-ray and 7-ray emission during only several hours, therefore 
this model requires quasi-continuous injection of energy into a single blob, or, 
most likely, a "multi-blob" scenario in which the observed radiation is a result 
of superposition of many short-live blobs continuously ejected from the central 
source. Then, the high state of the source like the extraordinary long outburst 
of Markarian 501 in 1997 could be associated with a dramatic increase of the 
rate of ejection/formation of 7-ray emitting blobs. 
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^ The minimum energy budget condition is achieved in the case of approximate 
equipartition between the protons and magnetic field, ~ W^. However, for 
Markarian 501 this condition would imply magnetic field of about 10 G, which only 
marginally could provide effective production of the proton-synchrotron radiation 
in the TeV regime. 
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